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ABSTRACT The Middle East respiratory syndrome (MERS) is proposed to be a zoonotic disease; however, the reservoir and mech- 
anism for transmission of the causative agent, the MERS coronavirus, are unknown. Dromedary camels have been implicated 
through reports that some victims have been exposed to camels, camels in areas where the disease has emerged have antibodies 
to the virus, and viral sequences have been recovered from camels in association with outbreaks of the disease among humans. 
Nonetheless, whether camels mediate transmission to humans is unresolved. Here we provide evidence from a geographic and 
temporal survey of camels in the Kingdom of Saudi Arabia that MERS coronaviruses have been circulating in camels since at 
least 1992, are distributed countrywide, and can be phylogenetically classified into clades that correlate with outbreaks of the 
disease among humans. We found no evidence of infection in domestic sheep or domestic goats. 

IMPORTANCE This study was undertaken to determine the historical and current prevalence of Middle East respiratory syndrome 
(MERS) coronavirus infection in dromedary camels and other livestock in the Kingdom of Saudi Arabia, where the index case 
and the majority of cases of MERS have been reported. 
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One hundred eighty laboratory-confirmed cases of human in- 
fection with Middle East respiratory syndrome coronavirus 
(MERS-CoV), 77 of them fatal, have been reported through 
30 January 2014 (1) following the identification of the index case 
in the Kingdom of Saudi Arabia (KSA) in September 2012 (2). The 
majority of infections have been identified in the KSA with lower 
numbers in Jordan, Qatar, Tunisia, and the United Arab Emirates. 
Although cases have also been reported in France, Germany, Italy, 
and the United Kingdom, all have been linked to the Middle East 
either by travel of the individuals infected through an area where 
MERS-CoV has been reported or by direct or indirect contact with 
others who have a travel history consistent with exposure in the 
Middle East (3). 

Clusters of human infection indicate that human-to-human 
MERS-CoV transmission can occur (4, 5). However, the origin of 
the infection in most cases remains unknown. Analysis of human 
MERS-CoV sequences by Cotten et al. has revealed the presence of 
at least three circulating genotypes within the KSA alone (6). Phy- 
logenetic analyses of 13 complete and 8 partial genome sequences 
enabled estimates of the timing and geographic origins of individ- 
ual viral clades. The authors proposed that MERS-CoV emerged 
in humans in 2011 and noted that sequence divergence between 



clades is consistent with several sporadic introductions of the virus 
into the human population, presumably from an animal reservoir. 

Efforts to identify an animal reservoir have focused on bats and 
camels. Bats harbor a wide range of betacoronaviruses (7); fur- 
thermore, bat cell lines display the MERS-CoV receptor, dipepti- 
dyl peptidase 4 (8), and can be experimentally infected. A short 
sequence fragment consistent with MERS-CoV was reported in a 
bat in Bisha, KSA, collected in close proximity to the home and 
workplace of the 2012 index case patient from whom the initial 
virus isolate was obtained (9). That same patient owned four pet 
dromedary camels (DC) . Serological analysis of those DC revealed 
the presence of antibodies reactive with MERS-CoV; however, no 
MERS-CoV sequences were found by PCR analysis of nasal or 
rectal swabs or serum. Additional human cases have been associ- 
ated with exposure to DC, and in some instances, investigators 
have described both serologic and genetic evidence of MERS-CoV 
infection in DC. Memish and coworkers reported PCR detection 
of MERS-CoV sequences in a DC with respiratory illness owned 
by an individual with MERS-CoV who had no history of contact 
with other infected humans (10). Haagmans et al. investigated an 
outbreak of the disease among humans on a Qatari farm and 
found MERS-CoV sequences in nasal swabs from 6 of 14 seropos- 
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TABLE 1 Samples collected in 2013 by animal species, geographic 
location, age group, and specimen type 



Animal species 


Location(s) 


Age group" 


No. 


Specimens'" 


DC 


Hofuf 


Juvenile 


19 


S, - 


B, N, R 


DC 


Hofuf 


Adult 


21 


S, 


B, N, R 


DC 




Tii vpn ilp 


21 


s, 


B, N, R 


DC 


Gizan 


Adult 


19 


s, 


B, N, R 


DC 


Taif 


Juvenile 


22 


S, - 


B, N, R 


DC 


Taif 


Adult 


19 


S, - 


B, N, R 


DC 


Tabuk 


Juvenile 


24 


s, 


8, N, R 


DC 


Tabuk 


Adult 


16 


s, 


B, N, R 


DC 


Riyadh 


Juvenile 


12 


c 

■3, - 


R NT D 
D, In, K 


DC 


Riyadh 


Adult 


8 


s, : 


3, N,R 


DC 


Unizah 


Juvenile 


6 


s, ; 


3, N,R 


DC 


Unizah 


Adult 


16 


s,: 


3, N,R 


Goat 


Unizah 


3 mo-2 yr 


31 


s,: 


3, N,R 


Goat 


Riyadh 


Unknown 


5 


s, . 


3, N,R 


Sheep, Barbari' 


Unizah 


3 mo-2 yr 


29 


s, . 


3, N,R 


Sheep, Harri 


Unizah, Riyadh 


3 mo-2 yr 


10 


s, . 


3, N,R 


Sheep, Najdi 


Unizah 


3 mo-2 yr 


21 


s,: 


B, N,R 


Sheep, Naimi 


Unizah 


3 mo-2 yr 


21 


s,: 


3, N,R 


Sheep, Sawakni 


Unizah 


3 mo-2 yr 


31 


s,: 


B, N,R 



" Juvenile animals were denned as being ^2 years of age; adults were denned as being 
>2 years of age. 

b Specimens collected from each animal: S, serum; B, blood; N, nasal swab; R, rectal 
swab. 

c Domestic sheep were separated into the breeds commonly found in the KSA: Barbari, 
Harri, Najdi, Naimi, and Sawakni. 

itive DC. Analysis of open reading frame la (ORFla) and frag- 
ments representing ORFlb, spike, and ORF4b revealed similarity 
but not identity to sequences obtained from the MERS-CoV- 
infected humans at the same farm. The authors provide evidence 
that MERS-CoV can infect DC but cautiously conclude that data 
are insufficient to determine whether the infection spread from 
DC to humans, from humans to DC, or via another host to both 
species (11). 

Several groups have reported serological reactivity with MERS- 
CoV or a closely related virus in DC in the Middle East (12-15). 
Reusken et al. found antibodies in 100% of50OmaniDCand 14% 
of 105 Canary Island DC but no seropositive northern European 
DC, domestic sheep, domestic goats, or domestic cattle (13). In 
two regions of the KSA, Hemida and colleagues detected antibod- 
ies to MERS-CoV in 90% of 310 DC but not in sheep, goats, cattle, 
or chickens. The seroprevalence was lower in DC < 1 year of age 
(72% versus 95%), suggesting widespread infection in early life 
(15). 

To determine the prevalence of MERS-CoV infection in DC 
throughout the KSA, we undertook a nationwide survey by using 
both serological and molecular methods. 

RESULTS 

Serum, whole blood, and rectal and nasal swabs were freshly col- 
lected from DC, sheep, and goats in November and December of 
2013 in the southwestern (Gizan), western (Taif), northwestern 
(Tabuk), eastern (Hofuf), and central (Unizah, Riyadh) regions of 
the KSA (Table 1). We also collected archived serum samples ob- 
tained from DC in 1992 through 2010 (Table2). Sera were initially 
tested for the presence of antibodies reactive with MERS-CoV by 
using a cell enzyme-linked immunosorbent assay (ELISA) based 
on Vero cells infected with MERS-CoV. Subsets of sera positive by 
ELISA were tested in Western blot assays that used extracts of Vero 



TABLE 2 Analysis of archived DC sera from the KSA from 1992 to 
2010 



Yr 


Location 


Age 
group 


No. 


% Seropositive 
(no. positive/total) 


1992 


Riyadh 


Adult 


1 


100 (1/1) 


1993 


Riyadh 


Adult 


2 


100 (2/2) 


1994 


Empty quarter 


Adult 


123 


93 (114/123) 


1996 


Riyadh 


Adult 


6 


100 (6/6) 


2004 


Riyadh 


Adult 


6 


100 (6/6) 


2009 


Riyadh 


Juvenile 


56 


72 (40/56) 


2009 


Rumah 


Adult 


26 


92 (24/26) 


2010 


Riyadh 


Juvenile 


21 


76 (16/21) 


2010 


Kharj 


Adult 


23 


91 (21/23) 



cells infected with MERS-CoV and a luciferase immunoprecipita- 
tion system (LIPS) assay based on recombinant MERS-CoV nu- 
cleoprotein. Potential serologic cross-reactivity with bovine coro- 
navirus (Bo-CoV) was addressed by testing for reactivity with Bo- 
CoV nucleocapsid protein by LIPS assay. The presence of viral 
nucleic acids in rectal and nasal swabs and a subset of serum and 
whole blood samples was assayed by reverse transcription- 
quantitative PCR (RT-qPCR) with primers targeting the upE and 
ORFla genome regions of MERS-CoV (16, 17). 

One hundred fifty (74%) of 203 DC sampled countrywide in 
2013 were found to have antibodies to MERS-CoV by ELISA. The 
prevalence of seropositivity was higher in adult DC (>2 years of 
age; 93/98, 95%) than in juvenile DC (<2 years of age; 57/104, 
55%) (P < 0.0001, x 2 test). The lowest prevalence of seropositive 
juveniles was found in the southwestern region, in proximity to 
the city of Gizan (1/21, 5%) (Fig. 1A). A higher prevalence of 
antibodies to MERS-CoV in older animals was also seen in sam- 
ples obtained in years prior to 2012. In 2010, the seroprevalence 
was 76% (16/21) in juveniles versus 91% (21/23) in adults in the 
central region (Riyadh, Kharj). In 2009, the seroprevalence was 
72% (40/56) in juveniles versus 92% (24/26) in adults (Riyadh, 
Rumah) (Table 2). Seroreactivity was also found in samples col- 
lected from DC as early as the 1990s: 1/1 (100%) in 1992, 2/2 
(100%) in 1993, 114/123 (93%) in 1994, 6/6 (100%) in 1996, and 
6/6 (100%) in 2003. Analysis of selected seropositive DC sera by 
Western immunoblot assay indicated two distinct reaction pat- 
terns that showed reactivity either to MERS-CoV spike glycopro- 
tein alone or to spike glycoprotein and nucleocapsid protein. 
These results were concordant with results of MERS-CoV nucleo- 
capsid LIPS assays (Fig. 2; see Table SI in the supplemental mate- 
rial). 

Potential serologic cross-reactivity to Bo-CoV was addressed 
by analyzing Bo-CoV nucleocapsid reactivity, which has been 
shown to be subgroup specific in CoVs (18). Overall, 17% (35/ 
203) of DC were positive for Bo-CoV, ranging from 3% in the 
southwest (Gizan) to 25% in the east (Hofuf) and 20% in adult 
animals versus 14% in juvenile animals; 2 animals (juveniles in 
Taif) were seropositive for Bo-CoV exclusively, while the remain- 
ing 16% (33/203) were reactive to Bo-CoV and MERS-CoV, and 
58% (1 17/203) were reactive to MERS-CoV alone (see Table SI in 
the supplemental material). In 2010, 47% (22/47) were Bo-CoV 
positive (57% [13/23] of adults and 38% [9/24] of juveniles, with 
1 adult exclusively positive for Bo-CoV), and in 2009, 20% were 
Bo-CoV positive (16/80; 8% [2/26] adults and 26% [14/54] juve- 
nile, with 2 juveniles exclusively positive for Bo-CoV). 
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Rectal and nasal swabs collected in parallel with serum samples 
from the same animals were assayed for MERS-CoV nucleic acids 
by RT-qPCR. Nucleic acids were most frequently detected in nasal 
swabs; rectal swabs were found positive in only three cases; in two 
of them, the nasal sample was also positive. The regional distribu- 
tion of PCR-positive animals is shown in Fig. IB. In contrast to 
serology, where MERS-CoV-reactive antibodies were more prev- 
alent in adults than in juveniles (95% versus 55%), MERS-CoV 
nucleic acids were found more frequently in juveniles (36/104, 
35%) than in adults (15/98, 15%) (P = 0.003, x 2 test). The five 
samples with >10 6 copies were all from juveniles, four of them 
from seronegative animals. The prevalence of PCR-positive DC 
ranged from 66% in Taif in the west to 0% in Gizan in the south- 
west. PCR analysis of a random selection of serum and whole 
blood samples collected from nasal or rectal swab PCR-positive, 
seropositive, and seronegative DC revealed no evidence of viremia 
(see Table SI in the supplemental material). These included 13 
adults and 29 juveniles phlebotomized in 2009, 15 adults and 14 
juveniles phlebotomized in 2010, and 8 adults and 13 juveniles 




H-9 

juvenile 
ELISA 
LIPSN»»* 
PCR neg 



H10 

juvenile 
ELISA 
LIPSN- 
PCR neg 



Tf-1 
adult 
EUSA 
LIPS N 
PCR ♦ nasal 



Tf-2 
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ELISA 
LIPS N 

PCR ♦♦» rectal 
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Rabbit anti- 

MERSCoV 
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Ctrl 

Rabbit anti- 
inactivated 
MERSCoV 



FIG 2 MERS-CoV proteins detected by Western blotting in sera from KSA 
DC. S, spike glycoprotein; N, nucleocapsid protein. 



phlebotomized in 2013. These animals included the five juveniles 
with the highest viral genome sequence load in nasal swabs. 

Serum samples collected from goats (« = 36) and sheep (n = 
112) in 2013 in the central region (Unizah, Riyadh) were not im- 
munoreactive with MERS-CoV but were immunoreactive with 
Bo-CoV(25%ofgoats, n = 36; 54% of sheep, « = 24). Nasal swabs 
from 36 goats and 78 sheep were negative in RT-qPCR assays for 
MERS-CoV upE. 

To test the validity of RT-qPCR results and determine phylo- 
genetic relationships of viral sequences found in DC in the KSA to 
previously reported sequences, we amplified and sequenced lon- 
ger regions of the spike, ORFlab, and nucleocapsid genes from 
RT-qPCR-positive samples (for the sequences of the primers used, 
see Table S2 in the supplemental material). Eleven of 13 swab 
samples with >10 5 copies in upE RT-qPCR yielded products for 
sequencing. No suitable products were obtained from samples 
with lighter viral sequence loads (<10 4 copies) (seeTableSl in the 
supplemental material). Phylogenetic analysis of a 1,044- 
nucleotide (nt) region of the spike gene and a 2,004-nt region of 
the ORF 1 ab gene indicated < 1 % divergence from previously pub- 
lished MERS-CoV sequences (Fig. 3) (GenBank accession no. 
KJ396756 to KJ396771). The nucleocapsid sequence was identical 
to the previously reported MERS-CoV sequences. 

DISCUSSION 

MERS-CoV is posited to be a zoonosis. However, the evolutionary 
history of MERS-CoV and the reservoirs and vectors for human 
infection remain obscure. Early anecdotal reports that some 
MERS-CoV victims had exposure to DC led to serologic investi- 
gation of DC in Spain and Oman (14), lordan (12), Egypt (13), 
and the KSA (15) that revealed antibodies to MERS-CoV. Defin- 
itive evidence that DC can be infected with MERS-CoV was ob- 
tained when viral sequences were detected in nasal swabs from DC 
sampled in close proximity to outbreaks of the disease among 
humans in Qatar (11) and Jeddah, KSA (10). Nonetheless, as 
noted by Nishiura and colleagues, current data do not fulfill the 
two criteria required to implicate DC as a significant reservoir 
species in the epidemiology of MERS-CoV (19), i.e., (i) that DC 
are sufficient to maintain MERS-CoV and (ii) that the presence of 
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FIG 3 Relationship of partial spike (A, 1,044 nt) and ORFlab (B, 2,009 nt) 
sequences derived from DC with selected human-derived MERS-CoV se- 
quences. Maximum-likelihood trees were inferred on the basis of p distance. 
The scale bars indicate the base difference per site. GenBank reference se- 
quence accession numbers are shown. 

DC is essential to the continuous transmission of infection. Re- 
sults presented here do not establish the latter; however, they do 
provide evidence for the former. 

Our study is the first comprehensive countrywide survey of DC 
from the KSA, the country with the most recorded MERS-CoV 
cases, to use both serological and molecular diagnostic methods. 
Analysis of specimens from western regions of the country, from 
Tabuk in the northwest, Taif in the west, and Gizan in the south- 
west, revealed regional differences. Although the seroprevalence 
was high in adults throughout the country at >80%, in juveniles, 
it ranged from 90% in the east to 5% in the southwest. The sero- 



prevalence in DC ^2 years of age was lower than that in older 
animals, confirming the results of Hemida et al. (15). Molecular 
analysis of nasal and rectal swab specimens indicated the highest 
prevalence of MERS-CoV sequences in DC in the west and north- 
west. Nasal swabs with heavy sequence loads (>10 5 copies) also 
clustered in the Taif region. A second sample collection in the west 
(Taif) separated from the first by an interval of 2 months con- 
firmed the presence of heavy viral sequence loads in nasal swabs 
collected from juvenile animals sampled in this area (data not 
shown). These findings suggest that continuous, longer-term sur- 
veillance is necessary to determine the dynamics of virus circula- 
tion in DC populations. Lower prevalence rates of both MERS- 
CoV and Bo-CoV were evident in samples from the southwest. 
This may relate in part to the enforcement of restrictions of live- 
stock movement in and out of Gizan Province implemented after 
the Rift Valley fever outbreak in 2000 but also to the generally 
lower DC population density in this region than in other regions 
of the KSA. 

Viral nucleic acids were more commonly detected in nasal 
swabs than in rectal specimens and were more frequent in juvenile 
than in adult animals. These findings, together with the absence of 
viremia and the known respiratory tract tropism of several other 
coronaviruses, suggest that airborne transmission is the most 
likely mode of MERS-CoV transmission. Although nucleic acid 
copy numbers were commonly highest in juvenile animals that 
were seronegative or had low antibody titers, positive findings 
were also obtained with specimens from highly seropositive and 
adult animals. 

Our findings in archived DC specimens, although restricted 
to serology, strongly suggest that MERS-CoV or a closely re- 
lated virus has been circulating in DC in the KSA for at least 2 
decades. Complete genomic sequences of MERS-CoV found in 
contemporary DC in the KSA are identical to sequences of 
viruses recovered from human MERS-CoV victims (unpub- 
lished data). Although we speculate that DC are potential res- 
ervoirs for human transmission, we cannot prove this relation- 
ship from the current data. Rigorous epidemiological 
investigation of the potential for exposure to DC in sporadic 
cases of MERS-CoV (those where there is no opportunity for 
human-to-human transmission) is required to test this model. 
If DC can be implicated, other questions will arise. Did MERS- 
CoV truly emerge as a human pathogen in 2012, or were cases 
of cryptic infection not appreciated because of a lack of suitable 
diagnostic tests? We may be able to address this conundrum by 
using archived human materials. If evidence of human MERS- 
CoV infections cannot be detected prior to 2012, we must en- 
tertain the possibility that mutation facilitated cross-species 
transmission. However, we see no path to address this possibil- 
ity absent access to historical DC respiratory tract specimens. 
The only archived DC specimens we have been able to locate 
are DC sera; our efforts to recover MERS-CoV sequences from 
camel blood have been unsuccessful. What are the roles of bats, 
if any, as reservoirs of MERS-CoV? These limitations notwith- 
standing, the most urgent public health concern, raised in work 
we and others have reported that focuses on DC infection, is to 
determine the role of these animals in sporadic human infec- 
tion. The evidence is clearly sufficient to support targeted in- 
vestigation of direct or indirect exposure to DC in the disease 
among humans. 
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MATERIALS AND METHODS 

Sample collection. Samples included DC, sheep, and goat sera; whole 
blood and nasal and rectal swabs freshly collected in 2013; and archived 
serum samples from 1992, 1993, 1994, 1996, 2004, 2009, and 2010. Two 
rectal and two nasal swabs were obtained from each animal. One rectal 
swab and one nasal swab were placed into RNAlater (Life Technologies, 
Carlsbad, CA), and one rectal swab and one nasal swab were placed into 
viral transport medium (Becton Dickinson, Franklin Lakes, NJ). All were 
stored at -80°C. 

ELISA. Vero (African green monkey kidney, ATCC CRL-1586) cells 
were maintained at the Integrated Research Facility (Frederick, MD) in 
Dulbecco's modified Eagle's medium (Corning Inc., Corning, NY) and 
10% fetal bovine serum. Cells were plated at a concentration of 4 X 10 4 / 
well in 96-well plates (catalog no. 3603; Corning). When cells were at or 
near confluence, they were infected with the Jordan strain of MERS-CoV 
(GenBank accession no. KC776174, MERS-CoV Hu/Jordan-N3/ 
2012[41]), kindly provided by Kanta Subbarao (National Institutes of 
Health, Bethesda, MD) and Gabriel Defang (Naval Medical Research 
Unit-3, Cairo, Egypt) at a multiplicity of infection of 1.0. At 24 h postin- 
fection, cells were fixed in 10% neutral buffered formalin or 4% parafor- 
maldehyde solution. After 24 h in fixative, plates were rinsed three times 
with phosphate-buffered saline (PBS) and placed in PBS for storage at 
4°C. Plates were loaded with infected and noninfected cells in alternating 
rows to generate a differential reading for each serum tested. Positivity was 
defined as an infected-cell optical density of >0.6 and >3X the 
noninfected-cell optical density. Test sera were diluted 1:3,000 in PBS- 
0.05% Tween 20-1% bovine serum albumin; secondary antibodies were 
rabbit anti-goat IgG (H+L)-horseradish peroxidase conjugate (1:3,000; 
Bio-Rad, Hercules, CA), rabbit anti-sheep IgG (H+L)-horseradish per- 
oxidase conjugate (1:3,000; Bio-Rad), and anti-llama IgG-horseradish 
peroxidase conjugate (1:10,000; Bethyl Laboratories, Montgomery, TX). 

Western blot. Extracts of noninfected Vero cells or Vero cells infected 
with MERS-CoV strain EMC were generated at Rocky Mountain Labora- 
tories, loaded onto discontinuous 3 and 7.5% SDS gels (Bio-Rad), and 
transferred onto nitrocellulose with iBlot Transfer Stacks (Invitrogen 
iBlot; Life Technologies). Lanes were loaded with alternating infected and 
noninfected extract samples, and a pair were cut for incubation with DC 
sera ( 1 :800 in blocking solution) after blocking of the membrane in PBS- 
0.05%Tween20-5%drymilkblocking solutionfor 1 h. Membranes were 
washed three times with PBS- 0.05% Tween 20 after a 2-h incubation with 
serum and then incubated for another 1.5 h with secondary antibody 
( 1 :7,000 in blocking solution; anti-llama IgG-horseradish peroxidase con- 
jugate; Bethyl Laboratories). Following three more washes, the mem- 
branes were developed with WesternSure premium chemiluminescent 
substrate (LI-COR, Lincoln, NE) and read on a C-DiGit Blot Scanner 
(LI-COR). 

LIPS assay. The nucleocapsid proteins of Bo-CoV and MERS-CoV 
were PCR amplified with primers introducing appropriate restriction 
sites for cloning into vector pREN-2 fused to the C terminus of the Renilh 
luciferase reporter (20). Sequence-confirmed construct DNA was purified 
from Escherichia coli cultures (Qiagen, Hilden, Germany), and transfected 
into COS-1 cells (African green monkey kidney, ATCC CRL-1650; 1 /xg; 
Lipofectamine; Invitrogen, Life Technologies). Cells were harvested at 
48 h posttransfection in lysis buffer (50 mM Tris [pH 7.5] , 100 mM NaCl, 
5 mM MgCl 2 , 1% Triton X-100, 50% glycerol, protease inhibitors), and 
the protein extract was clarified by two rounds of centrifugation at 12,000 
X g. The target concentration was determined in relative light units 
(RLU), and approximately 10 X 10 6 RLU were incubated with test serum 
( 1: 100) in a final volume of 100 iaI buffer A (50 mM Tris [pH 7.7] , 100 mM 
NaCl, 5 mM MgCl 2 , 1% Triton X-100) per well in 96-well plates (catalog 
no. 249944; Thermo/Nunc, Waltham, MA). After 1 h of incubation at 
room temperature, protein A/G beads (catalog no. 53 135; Pierce, Junction 
City, OR) were added and the mixtures were transferred to 96-well filter 
plates (catalog no. MSBVN1B50; Millipore, Billerica, MA) for another 
hour of incubation at room temperature with shaking. Bead-bound anti- 



gen was washed eight times with buffer A and three times with PBS (Tecan 
Hydroflex, Maennedorf, Switzerland) and then read with coelenterazine 
substrate (Renilla luciferase assay system; Promega, Madison, WI) on a 
Centro LB960 luminometer (Berthold, Bad Wildbad, Germany). 

Nucleic acid extraction and PCR. Total nucleic acids were extracted 
from nasal swabs, serum, and whole blood on a QiaCube with Cador 
Reagent kits (Qiagen) or RNeasy Reagent kits for extraction of RNA from 
rectal swabs. Real-time qPCR used a OneStep Real-Time qPCR buffer 
(Invitrogen, Life Technologies) and primer/probes upE and ORFla (16, 
17). Products for sequencing were generated by RT-PCR. cDNA was re- 
verse transcribed with Superscript III and random hexamer primers. PCR 
was performed with Amplitaq Gold (Life Technologies) and primers de- 
signed to amplify a 1,044-nt region of the spike gene (heminestedPCR), a 
913-nt region of the N gene (nested PCR) or a 2,004-nt region of the 
ORFlab region (heminested PCR). For the sequences of the primers used, 
see Table S2 in the supplemental material. Products were purified by 
agarose gel electrophoresis and with QIAquick Gel Extraction kits (Qia- 
gen) and subsequently sequenced on both strands by the dideoxynucle- 
otide chain termination method (GeneWiz, South Plainfield, NJ). 

SUPPLEMENTAL MATERIAL 

Supplemental material for this article may be found at http://mbio.asm.org 
/lookup/suppl/doi:10.1128/mBio.00884-14/-/DCSupplemental. 

Table SI, DOCX file, 0.1 MB. 

Table S2, DOCX file, 0.1 MB. 
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